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Quantum-mechanically correlated (entangled) states of many particles are of interest in quan¬ 
tum information, quantum computing and quantum metrology. Metrologically useful entangled 
states of large atomic ensembles have been experimentally realized p-10, but these states display 
Gaussian spin distribution functions with a non-negative Wigner function. Non-Gaussian entangled 
states have been produced in small ensembles of ionsfTTl [12] , and very recently in large atomic 
ensembles [I3hl 5) . Here, we generate entanglement in a large atomic ensemble via the interaction 
with a very weak laser pulse; remarkably, the detection of a single photon prepares several thou¬ 
sand atoms in an entangled state. We reconstruct a negative-valued Wigner function, an important 
hallmark of nonclassicality, and verify an entanglement depth (minimum number of mutually en¬ 
tangled atoms) of 2910 d= 190 out of 3100 atoms. This is the first time a negative Wigner function 
or the mutual entanglement of virtually all atoms have been attained in an ensemble containing 
more than a few particles. While the achieved purity of the state is slightly below the threshold 
for entanglement-induced metrological gain, further technical improvement should allow the gener¬ 
ation of states that surpass this threshold, and of more complex Schrodinger cat states for quantum 
metrology and information processing. More generally, our results demonstrate the power of her¬ 
alded methods for entanglement generation, and illustrate how the information contained in a single 
photon can drastically alter the quantum state of a large system. 


Entanglement is now recognized as a resource for se¬ 
cure communication, quantum information processing, 
and precision measurements. An important goal is the 
creation of entangled states of many-particle systems 
while retaining the ability to characterize the quantum 
state and validate entanglement. Entanglement can be 
verified in a variety of ways, with one of the strictest 
criteria being a negative-valued Wigner function |T6l ITT], 
that necessarily implies that the entangled state has a 
non-Gaussian wavefunction. To date, the metrologically 
useful spin-squeezed states [TUTO] have been produced in 
large ensembles. These states have Gaussian spin distri¬ 
butions and therefore can largely be modeled as systems 
with a classical source of spin noise, where quantum me¬ 
chanics enters only to set the amount of Gaussian noise. 
Non-Gaussian states with a negative Wigner function, 
however, are manifestly non-classical, since the Wigner 
function as a quasiprobability function must remain non¬ 
negative in the classical realm. While prior to this work 
a negative Wigner function had not been attained for 
atomic ensembles, in the optical domain, a negative¬ 
valued Wigner function has very recently been measured 
for states with up to 110 microwave photons [ 18 k Another 
entanglement measure is the entanglement depth [19], i.e. 
the minimum number of atoms that are demonstrably, 
but possibly weakly, entangled with one another. This 
parameter quantifies how widely shared among the par¬ 
ticles an entangled state is. For a state of an ensemble 
characterized by collective measurements, the entangle¬ 
ment depth depends sensitively on the proximity of the 
state to the ideal symmetric subspace of all particles. 


The largest entanglement depth verified previously has 
been 170 out of 2300 atoms for a spin-squeezed state [6]. 
and very recently 13 out of 41 atoms for a non-Gaussian 
state [13]. 

Here we generate entanglement in a large atomic en¬ 
semble by detecting a single photon that has interacted 
with the ensemble [20]. An incident vertically polarized 
photon experiences a weak random polarization rota¬ 
tion associated with the quantum noise of the collective 
atomic spin. The detection of a horizontally polarized 
emerging photon then heralds a non-Gaussian entangled 
state of collective atomic spin (Fig. 1) with a negative¬ 
valued Wigner function of —0.36 =b 0.08, and an entan¬ 
glement depth of 90% of our ensemble containing several 
thousand atoms. 

The pertinent atom-light interaction is enhanced by 
an optical cavity, into which we load N a = 3100 db 300 
laser-cooled 87 Rb atoms (Fig. la). The atoms are pre¬ 
pared in the 5 Si/ 2,E = 1 hyperfine manifold, such that 
each atom i can be associated with a spin f i: and the 
ensemble with a collective-spin vector S = Yhi fi- After 
polarizing the ensemble (S z « S) by optical pumping, the 
collective spin state is rotated onto the x axis by means 
of a radiofrequency tt/2 pulse. This (unentangled) ini¬ 
tial state that is centered about S z = 0 with a variance 
(A S z ) 2 = S/2 is known as a coherent spin state (CSS). 
In our experiment, the atoms are non-uniformly coupled 
to the optical mode used for state preparation and de¬ 
tection, but the relevant concepts can be generalized to 
this situation, as discussed in Methods. 

Probe light resonant with a cavity mode and detuned 
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FIG. 1: Scheme for heralded entanglement generation in a large atomic ensemble by single-photon detection, (a) Incident 
vertically polarized light experiences weak polarization rotation due to atomic quantum noise, and the detection of a hori¬ 
zontally polarized transmitted photon heralds an entangled state of collective atomic spin. An optical resonator enhances the 
polarization rotation and the heralding probability, (b) Atoms in the 551/2, F — 1 hyperfine manifold are coupled to the excited 
5 P 3/2 manifold via linearly polarized light, decomposed into two circular polarization components |cr^) that interact with the 
atomic ground-state populations. The outgoing polarization state of the light reflects the quantum fluctuations between the 
\5S 1/2 F = l,ra = ±1) magnetic sublevels. 


from the 87 Rb D 2 transition is polarization analyzed 
upon transmission through the cavity. The vertical po¬ 
larization state of each photon in the incident laser pulse 
|v) = (|cr + ) + \<j~))/\/2 can be decomposed into two cir¬ 
cular polarization components Icr^) that produce oppo¬ 
site differential light shifts between the atomic magnetic 
sublevels \m = ±1). Hence a |cr ± ) photon causes a pre¬ 
cession of the collective spin vector S in the xy plane 
by a small angle ±0 (see Methods), and we denote the 
corresponding slightly displaced CSS by |d=0). Then the 
combined state of the atom-light system after the passage 
of one photon can be written as [20] 

\i>) oc\a + )\+(j>)+ \a~)\-(j>). ( 1 ) 

Conversely, atoms in the states \m = ±1) cause different 
phase shifts on the cr ± photons, resulting in a net rotation 
of the photon linear polarization if the states |m = ±1) 
are not equally populated. Then the atomic quantum 
fluctuations between \m = ±1) in the CSS randomly ro¬ 
tate the polarization of the input photons |v), giving rise 
to a nonzero probability oc 0 2 for an incident |v) photon 
to emerge in the polarization |h) = (|cr + ) — |cr - ))/\/2, 
orthogonal to its input polarization. The detection of 
such a “heralding” photon projects the atomic state onto 
(h^) oc |0) — |—0), which is not a CSS, but an entangled 
state of collective spin, namely, the first excited Dicke 
state [21] I'lpi) along x (Fig. la). In contrast, if the 
photon is detected in its original polarization |v), the 
atomic state is projected onto (v^) oc |0) + |—0), a state 
slightly spin squeezed [li and essentially identical to the 


input CSS. Thus the entangled atomic state |t/q) is post- 
selected by the detection of the heralding photon |h). 

From a different perspective, the entangled state is gen¬ 
erated by a single-photon measurement event. The inci¬ 
dent photon undergoes Faraday rotation by an angle 7 ? 
proportional to the collective spin along the cavity axis, 
S z , that exhibits quantum fluctuations around (S z ) = 0. 
Since detection of the outgoing photon in |h) is only pos¬ 
sible if S z % 0, such detection excludes values of S z near 
0 from the spin distribution [20], and biases the collective 
spin towards larger values of \S Z \. This creates a “hole” 
in the atomic distribution near S z = 0, as seen in Fig. 
la. 

The mean photon number in the incident laser pulse 
k ~ 210 is chosen such that the probability for one pho¬ 
ton to emerge in heralding polarization |h) is p ~ 0.05 <C 
1. This ensures a very small probability oc p 2 for pro¬ 
ducing a different entangled state | 02 ) heralded by two 
photons [20], a state which, due to our photon detection 
efficiency of q = 0.3 < 1, we would (mostly) mistake 
for IVh). This admixture of |02) to the heralded state 
is suppressed by a factor of 3p(l — q) « 0.1. Further 
state imperfection arises from false heralding events due 
to residual polarization impurity of the probe beam (in¬ 
dependent of the atoms) of ~ 3 x 10 -5 = 0.1p/k, adding 
an admixture of about 10% of the CSS to the heralded 
state. 
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FIG. 2: Collective-spin distribution of atomic state heralded 
by one photon, (a-b) Measured photon distributions g(np) 
for no heralding photon detected (blue squares), and for one 
heralding photon detected (red circles), for rotation angles 
(a) /3 = 0, (b) /3 = 7r/4, (c) /? = 7r/2, (d) /? = 37r/4. In- 
set: Logarithmic representations of the same data. In the 
ideal case, the ratio for the heralded state and the CSS is 
(richer/(rip)css = (#|)her/(£|)css = 3 for any angle /3, and 
we measure (np}her/(np)css = {2.7 zb 0.2, 2.2 ± 0.2, 2.4 ± 
0.2, 2.1 =b 0.1} for /3 m. {0, \^}. For each /3, the blue 
and red data sets represent approximately 1.5 x 10 4 and 200 
experiments, respectively. The solid blue and the dashed red 
curves are predictions without any free parameters, calcu¬ 
lated from first principles and the separately measured atom 
number, for the CSS and the perfect first Dicke state, respec¬ 
tively. The solid red line corresponds to the simultaneous fit 
to all measurement angles /3, i.e. the reconstructed density 
matrix. Error bars indicate 1 standard deviation (s.d.) (e- 
h) Reconstructed collective spin distributions of the heralded 
state (red) for rotation angles (e) f3 — 0, (f) [3 = 7r/4, (g) 
/ 3 = 7 t/2 , (h) [3 = 37r/4. The spin distributions of the CSS 
(blue) are for reference. The horizontal axis S z is expressed in 
terms of the effective atom number[4] N = (2/3) N a = 2100, 
obtained by weighting each atom with its coupling strength to 
the standing-wave probe field inside the cavity, such that the 
experimentally measured spin fluctuation (AS'*) 2 of the CSS 
via its interaction with the probe light satisfies the standard 
relation (AS Z ) 2 = S/2 = NF/2 for spin F atoms (see Meth¬ 
ods). The shaded area indicates the statistical uncertainty of 
1 s.d. The spin distribution in Fig. 2f shows no “hole” in the 
middle due to lower quality of data for this measurement run 
13 = tt/4. 


In order to reconstruct the collective-spin state gener¬ 
ated by the heralding event, we rotate the atomic state 
after the heralding process by an angle (3 = 0, f, f, 
about the x axis before measuring S z . (Thus (3 = 0 
corresponds to measuring S z , j3 = 7t/ 2 corresponds to 
S y , etc.) The measurement is performed by applying a 
stronger light pulse in the same polarization-optimized 
setup used for heralding. As the Faraday rotation angle 
i3 <C 1 is proportional to S z , and the probability for de¬ 
tecting |h) photons is proportional to $ 2 , the measured 
probability distribution of |h) photon number, g(np ), re¬ 
flects the probability distribution of Fig. 2a-d show 
that a single heralding photon substantially changes the 
spin distribution towards larger values of (S|). We fur¬ 
ther verify that the heralded state remains (nearly) spin 
polarized with a contrast of C = 0.99^o‘.o2? the same as 
for the CSS within error bars (Fig. 3a). 

From the photon distributions g(np) we can recon¬ 
struct the density matrix p rnn in the Dicke state basis[21] 
along x, where \n = 0) denotes the CSS along x, | n = 1) 
the first Dicke state, \n = 2) the second Dicke state, etc. 
From the density matrix we obtain the Wigner function 
W{6 , 0) on the Bloch sphere[22] (Fig. ??). To accurately 
determine the Wigner function value on the axis, W(6 = 
^,0 = 0) = ^Z n {—l) n Pnm that depends only on the pop¬ 
ulation terms p nn , we average the photon distributions 
g{np) over four angles (3 and thereby reduce the fitting 
parameters to just p nn ,n < 4. This is equivalent to con¬ 
structing a rotationally symmetric Wigner function from 
the angle-averaged marginal distribution[17 . We obtain 
poo = 0.32±0.03, pi i = 0.66±0.04 with negligible higher- 
order population terms, giving FF(|, 0) = —0.36 ± 0.08, 
to be compared to FF(|, 0) = —1 for the perfect first 
Dicke state. 

We can also fit the density matrix including the coher¬ 
ence terms simultaneously to g(np) for all four angles /?, 
without angle-averaging. Since the photon distributions 
g(np) depend only on they determine only the even 
terms of the density matrix, i.e., p mn where ra+n is even, 
and contain no information about the odd terms. If we 
calculate W(‘- 1,0) from the density matrix without angle¬ 
averaging, we find FF^O) = —0.27 ± 0.08, within error 
bars consistent with the angle-averaged value. In order 
to display the Wigner function, we bound the odd terms 
(ra + n odd) by verifying that the heralding process does 
not displace the state relative to the CSS (see Methods). 
Therefore we set the odd terms to zero, and display the 
resulting density matrix and corresponding Wigner func¬ 
tion in Fig. 3b-d. The spin distributions f(Sp) obtained 
from this density matrix are shown in Fig. 2e-h. 

In order to quantify the minimum number of mutually 
entangled atoms, we use a criterion derived in Ref. m 
that establishes entanglement depth as a function of the 
populations poo an d pn- From this criterion, generalized 
to the case of non-uniform coupling to the measurement 
light field (see Methods), we deduce an average entangle- 
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FIG. 3: Reconstruction of the heralded many-atom entangled state, (a) Normalized spin component S z /S measured in a 
Ramsey sequence, as a function of the phase of the second Ramsey 7r/2 pulse, for the CSS (blue squares) and the heralded 
state (red circles). The fit (red line) shows a contrast of 0.99^°'^ for the heralded state, within error bars the same as the 
contrast 0.995 ± 0.004 of the CSS. The negligible contrast reduction is expected given that we send only 210 photons into the 
system at large detuning from atomic resonance, (b) Reconstructed Wigner function VF(0,0) for the heralded state on the 
Bloch sphere [22] with a radius given by the effective atom number N — 2100. 6 is the polar angle with respect to z and cj) 
is the azimuthal angle with respect to x. The first excited Dicke state and the CSS have VF(|, 0) = —1 and 0) = 1, 

respectively. To provide a reference scale for the size of the negative region, the black dashed line is the contour at which 
the CSS has a Wigner function value 1/e. (c),(d) Real and imaginary parts of the reconstructed density matrix elements, in 
the Dicke state basis along x , for the heralded state, (e) Entanglement depth criterion[13] for the heralded state, plotted in 
terms of density matrix elements poo and p\\. The red shaded region represents the 1 s.d. confidence region for the heralded 
state. Lines represent boundaries for /e-particle entanglement in terms of atom number N a ; a state with pn greater than such 
a boundary displays at least /c-particle entanglement. States falling within the blue shaded region are not provably entangled 
by the used criterion. The hatched area indicates the unphysical region where the density matrix trace would exceed unity. 


ment depth of J\f a = 2910 ± 190 out of N a = 3100 atoms 
(Fig. 3e) using the angle-averaged density matrix. Our 
results represent the first experimental verification of the 
mutual entanglement shared by virtually all atoms in an 
ensemble that contains more than a few particles. 

The above results demonstrate that even with lim¬ 
ited resources, i.e. weak atom-photon coupling, herald¬ 
ing schemes can be used to boost the effective inter¬ 
action strength by a large factor, enabling the pro¬ 
duction of highly entangled states [20] [23]. Further¬ 
more, by repeated trials and feedback the entangle¬ 
ment generation can be made quasi-deterministicj24j|25]. 
Our approach is related to other heralded schemes 
for quantum communication [24U27] and entangled-state 
preparation [28H3Q] . and it would be interesting to gen¬ 
eralize the present analysis to infer characteristics of the 
atomic state from the measured optical signals in those 
experiments. We note that the same first Dicke state was 
created in an ensemble of up to 41 atoms with a scheme 
that uses many heralding photons in a strongly coupled 
atom-cavity system[ 13 . In our system, the maximum 
atom number of ~ 3000 is set by the accuracy of the spin 
rotation, and can be increased by two orders of mag¬ 


nitude by better magnetic-field control [lOj. The state 
purity pn can probably be further improved by reduc¬ 
ing the heralding probability, and a value of pu > 0.73 
would be required for the Fisher information 14] to ex¬ 
ceed that of the CSS, and enable metrological gain of 
up to 3 dB. The detection of two or more photons pre¬ 
pares Schrodinger cat states [20] of the atomic ensemble 
with more metrological gain. We expect that heralded 
methods can generate a variety of nearly pure, complex, 
strongly entangled states that are not accessible by any 
other means at the present state of quantum technology. 
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PHOTON POLARIZATION ROTATION BY 
ATOMIC SPIN 


Probe laser light red-detuned by A 0 /( 27 r) = —200 MHz 
from the 87 Rb transition 5 2 S'i/ 2 , F = 1 to 5 2 P 3 / 2 , F' = 0 
is sent through an optical cavity containing the atomic 
ensemble. We first consider the case where all the atoms 
are coupled with equal strength to the probe light. For 
detuning A much larger than the excited state linewidth 
r/( 27 r) =6.1 MHz, the excited state manifold can be 
adiabatically eliminated. The vector component of the 
ac Stark shift is described by the Hamiltonian 


H 

h 



( 2 ) 


where J z = ^(a^a + — ala_), with a± the annihila¬ 
tion operators for photons with cr ± circular polariza¬ 
tions. Here 2 g is the effective single-photon Rabi fre¬ 
quency taking into account the multiple transitions from 
5 2 Si/2, F = 1 to 5 2 P 3/2 , F' = 0,1, 2, given by 


s 2 = (< 7?;?) 2 + (^ 


1,0x2 


2,0x2 


2,2x2 


( 3 ) 


where 2 is the single-photon Rabi frequency be¬ 
tween the ground state | F = 1, m) and the excited state 
| F',m'). As Ao is comparable to the hyperfine split¬ 
tings of the 5 2 P 3/2 excited states, the interaction strength 
g 2 / A is given by 
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A 0 A 0 — Ai A 0 — Ai — A 2 A 0 — Ai — A 2 


where Ai/(27r) = 72 MHz is the hyperfine splitting be¬ 
tween the F' = 0 and F' = 1 manifolds, A 2 /(27r) = 157 
MHz between F' = 1 and F' = 2, and A/(2n) = —150 
MHz is the effective detuning when Aq/(27t) = —200 
MHz. The value g 2 / A for our experiment is 2tt x 0.7 
kHz. 

This vector shift © gives rise to a J z -dependent Lar- 
mor precession of the atomic collective spin S in the xy 
plane. Consider one Icr^) photon passing through the 
optical cavity and causing the atomic spin to precess by 
phase ±0. The characteristic atom-photon interaction 
time is 2/ft, where ft is the cavity linewidth, therefore the 
atomic phase is given by[lj[2] (j) = g 2 /(An) = 7^ v r/(4A), 
where the cavity cooperativity g v = 4g 2 /(/T) = 0.07. 


Another way to think of the Hamiltonian © is that the 
atomic spin component S z causes different phase shifts 
on the photon cr + and a~ components, resulting in a 
rotation of the linear polarization of the light. The po¬ 
larization rotation angle $ = (g 2 /A) (S z / 2) (2/ft) = (j)S z . 

In general, the incident light can introduce Raman 
transitions between different magnetic levels in the F = 1 
ground state manifold. We apply a bias magnetic field of 
4.7 G along the cavity axis to introduce a Zeeman shift 
between the magnetic levels, so that the Raman coupling 
is off-resonant. The Larmor frequency is cjl/(27t) = 3.3 
MHz, larger than the cavity linewidth k/(2it) = 1.0 MHz, 
so that the Raman coupling can be neglected. There is 
also an unimportant scalar light shift, as well as a tensor 
light shift that gives rise to squeezing that is negligible 
for our experimental conditions. 


EXPERIMENTAL DETAILS 

We load an ensemble of 87 Rb atoms, cooled to T = 
50 fi K, into a medium-finesse optical cavity (cavity finesse 
T = 5600, linewidth k/{2tt) = 1.0 MHz, cooperativity 
770 =0.2 at an antinode on a transition with unity oscil¬ 
lator strength). The atoms are confined on the cavity axis 
by a far-detuned optical dipole trap at 852 nm with trap 
depth U/h = 20 MHz. Characteristics of the optical cav¬ 
ity at the 780 nm probe laser wavelength and the 852 nm 
trap laser wavelength are summarized in Extended Data 
Table 1. One Glan-Taylor polarizing beamsplitter (Thor- 
labs GT5) purifies the polarization of probe light entering 
the cavity, while a second polarizing beamsplitter after 
the cavity allows us to measure the rotation of the probe 
light due to the atomic projection noise. Two Single Pho¬ 
ton Counting Modules (SPCMs, models SPCM-AQRH- 
14-FC and SPCM-AQR-12-FC) are placed at the trans¬ 
mitting and reflecting ports of the polarizing beamsplit¬ 
ter to detect the photons. Due to the fiber coupling and 
finite SPCM detection efficiency at 780 nm, the overall 
quantum efficiency of the detection process is q = 0.3. 


DEFINITION OF EFFECTIVE ATOM NUMBER 

Atoms are optically confined at the antinodes of the 
852 nm trap laser standing wave. The 780 nm probe 
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Parameter 

A=780nm 

A=852nm 

Mirror separation 

L 

26.62(1) mm 

Mirror curvature radius 

R 

25.04(2) mm 

Free spectral range 

<*>fsr/ (27r) 

5632.0(2) MHz 

Transverse mode spacing 

6 > a /(2 n) 

226.3(3) MHz 

Linewidth 

Kj(2n) 

1.01(3) MHz 

135(2) kHz 

Finesse 

h 

5.6(2)x 10 3 

4.2(l)x 10 4 

Mode waist 

WA 

56.9(4) pm 

59.5(5) pm 

Antinode cooperativity 

Vo, A 

0.203(7) 

1.65(4) 


Extended Data Table I: Resonator parameters. The mode 
waists are calculated at the position of the atoms. Outside 
this table, all resonator values refer to the probe wavelength 
A — 780 nm. 


light in the cavity forms a standing wave that is incom¬ 
mensurate with the trap standing wave. Consequently, 
the atoms experience spatially varying couplings to the 
probe light and rotate the probe photon polarization by 
different amounts. For an atom at position 2 on the cav¬ 
ity axis, the cooperativity is rj(z) = r] v sin 2 (kz). When 
N a atoms are prepared in a CSS, the atomic projection 
noise gives rise to fluctuations of the photon polariza¬ 
tion rotation. The measured variance of the polarization 
rotation is proportional to ^{p 2 (z)) where averaging 
is performed over the position z. This variance differs 
by a factor of order unity from that of a CSS consist¬ 
ing of N a atoms uniformly coupled to the light. As de¬ 
scribed in a previous paper [31, we introduce the effective 
atom number N and the effective cavity cooperativity r] 
to satisfy two conditions: that the experimentally mea¬ 
sured variance equals that of N uniformly coupled atoms, 
^f{rj 2 } = y 7 ? 2 ; and that the total amount of interaction 
between the atomic ensemble and the probe light is the 
same, i.e., N a (r] z ) = Nrj. To satisfy these two conditions 
we define the effective atom number N = |7V a and the 
effective cavity cooperativity r] = ^rj v . This re-scaling al¬ 
lows direct comparison with the well-known expressions 
for the uniformly coupled CSS. 

As in the main paper and the rest of Methods, S z 
refers to the collective spin of an ensemble containing N 
effective atoms, and therefore the atomic spin precession 
phase for each transmitting cavity photon is given by 
(j) = rjT / (4A) = (3/4)?7 v r/(4A). Note that this value 77 = 
0.05 < 1 corresponds to the weak atom-cavity coupling 
regime. For our parameters, = 5 x 10 4 < <ficss = 
1.5 x 10 -2 where <ficss — 1/(25) is the angular rms 

width of the CSS. 


CHOICE OF THE HERALDING PHOTON 
NUMBER 

The heralding light must be weak enough that it 
does not introduce substantial decoherence of the de¬ 
sired atomic state. The fundamental shot noise between 
the cr + and a~ circular polarization components of the 
heralding light gives rise to phase broadening of the 
atomic state, which limits the purity of the heralded en¬ 
tangled state. To measure the phase broadening, herald¬ 
ing light pulses with variable photon number are sent 
into the cavity, and the variance A S' 2 is measured by ap¬ 
plying a radiofrequency n/2 pulse to rotate the atomic 
state about the x direction before measuring A S' 2 . Ex¬ 
tended Data Fig. 1 shows the measured atomic state 
variance A S 2 as a function of the photon number in the 
heralding light, in agreement with the predicted linear 
dependence. The heralding photon number is thus cho¬ 
sen to be ^ 210, with corresponding herald detection 
probability qp = 1.5%, to give fairly small phase broad¬ 
ening. Lower heralding photon number results in a purer 
heralded state, but at the expense of a lower heralding 
and state generation probability. 


heralding light input photon number 
0 500 1000 1500 2000 2500 



herald detection probability pq 


Extended Data Figure 1: The measured atomic state variance 
A Sy as a function of the heralding light photon number and 
corresponding probability qp of detecting one photon. The 
solid red line is the prediction for A Sy broadened by the pho¬ 
ton shot noise of the heralding light. The dashed black line 
shows the CSS variance for 2030 F — 1 effective atoms used 
in this measurement. 


RELATION BETWEEN THE SPIN 
DISTRIBUTION f(Sp) AND THE MEASURED 
PHOTON DISTRIBUTION g[np) 

To measure the atomic state spin distribution, mea¬ 
surement light with the same polarization |v) as the 
heralding light is sent through the atoms, and the number 
of photons with the orthogonal polarization |h) is mea- 
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sured. The measurement light contains a large number 
of input photons ri[ n = 1.7 x 10 4 to perform destruc¬ 
tive measurements with good signal-to-noise ratio. The 
photon polarization is rotated by a small angle P = <j)S z 
and the probability for each photon to emerge in |h) is 
'd 2 . For a given number of input photons n i n , the aver¬ 
age number of detected photons with |h) polarization is 
(n) = qni n ((/)S z ) 2 , where q is the overall quantum effi¬ 
ciency. Therefore, a spin distribution f(S z ) is mapped 
to a measured photon distribution g(n). For a given S z , 
the detected photons follow a Poisson distribution with 
the mean number (n), and the probability to measure 
exactly n photons is given by 

P(fl)S z ) — exp[ qTliny^Pz) ] ~y • (4) 

For an atomic state with the spin distribution f(S z ), the 
photon distribution g{n) is given by 

g(n) = J2f(S z )P(n,S z ) 

s z 

ST' f(C \ r (AC \2i [Q n in{^S z ) 2 ] n 

= 2^/(^) exp [ _ ^ nin v^) ]-j-• 

n • 

(5) 

In order to measure the spin along a general direction, 
the atomic spin is rotated by an angle /? with a radiofre¬ 
quency pulse prior to detection. Replacing S z by Sp in 
equation © we write the relation between the spin dis¬ 
tribution f(Sp) and the measured photon distribution 
g(n 0 ) as 

9( n p) = J2f(S 0 )P(n 0 ,S 0 ) 

S/3 

= '^2f(S 0 )exp[-qn in ((f>S 0 )‘ 2 ] ^ qrim ^ S ^ l —. 

Q n P m 

(6) 

CHOICE OF THE MEASUREMENT PHOTON 
NUMBER 

The measurement photon number is chosen to optimize 
the readout quality. Extended Data Fig. 2 illustrates the 
dependence of readout on the input measurement photon 
number ni n by showing how the reconstructed distribu¬ 
tions f(S z ) change as n- in is varied (the method of recon¬ 
struction is discussed later). When the photon number 
is small, there is large detection noise due to photon shot 
noise, reflected as the large error band. With increasing 
photon number, the photon scattering by atoms into free 
space increases and the atomic state is more strongly 
perturbed, therefore the “dip” at S z = 0 becomes less 
distinct. To balance these two competing effects, the op¬ 
timized atomic-state-measurement photon number is set 
to 1.7 x 10 4 . 


SUBTRACTING BACKGROUND PHOTON 
COUNTS 

Due to the residual polarization impurity of the mea¬ 
surement light, there are a small number of background 
photon counts even when there are no atoms. The back¬ 
ground counts account for about 4% of the photon sig¬ 
nal of the heralded state. We independently measure 
the background photon distribution and subtract it from 
the directly measured atomic signal to obtain g(np). If 
we were not to correct for these background counts, we 
would overestimate the density matrix population pn by 
10 %. 


RECONSTRUCTION OF THE DENSITY 
MATRIX 

Using the measured photon distributions g(np) for all 
four angles /3 = 0, tt/4, 7t/2, 37t/ 4, the density matrix p of 
the heralded state can be reconstructed. 

As the entangled state maintains 0.991 q!o 2 contrast, 
the length of the total spin S « N and we can express 
the density matrix in the basis of Dicke states | m) x along 
the x direction 

P = Poo |0)a; (Ola; + Pn| l)x(lU + POl|0)x(lU + PlO | l)cc(0|x 

+P22 |2)a;(2|cc + P02 |0) x {2\ x + P20 |2) x (0|^ +- (7) 

The spin distribution f(Sp) can be written as a func¬ 
tion of atom number N and the density matrix elements 
Poo>Pn> e t c: 

f{Sp,p,N) = (Sp\p\Se) 

= p 00 G(0, Sp)G*{ 0, Sp) + pnG(l, Sp)G*( 1, Sp) 
+ P oiG(0, Sp)G*( 1, Sp) + pi 0 G(l, Sp)G*( 0, Sp) 
+P 22 G& Sp)G*( 2, Sp) + po 2 G( 0, Sp)G*( 2, Sp) 
+P2oG(2, Sp)G*(0, Sp) + — (8) 

Here G(m,Sp ) = (Sp\m) x is the wavefunction of Dicke 
state | m) x in the representation of spin component Sp 
and is given by 

1 / 1 X 1 / 4 

G(m, S/ 3 , N) = - — x 

P y/2™m\ \ttNJ 

xe im 0 -sj/(2N) Hm , (9) 

where H m (x) is the rath order Hermite polynomial and 
N is the atom number. Using equation (J6|, we write the 
theoretically predicted photon distribution gth{np) as a 
function of the density matrix p, atom number N and 
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Extended Data Figure 2: Dependence of the reconstructed distribution of collective spin S z on the measurement photon number, 
as illustrated by reconstructed spin distributions for photon numbers (a) 0.5 x 10 4 , (b) 1.1 X 10 4 , (c) 1.7 x 10 4 , (d) 2.7 x 10 4 , 
(e) 3.6 X 10 4 . Blue lines correspond to the CSS and red lines correspond to the heralded states. The shaded area indicates an 
uncertainty of 1 standard deviation. 


input photon number n\ n 

9th(P' ft1 P 5 N\ ^in) ^ ^ p , N^Pjjlp^ Sp') 

= ^2 fth{S/3, p, N) exp[-qn in (Sp(p) 2 ] \ —. 

a, n ?' 

( 10 ) 

We independently measure the input photon number ri[ n 
and find the atom number N by fitting the photon dis¬ 
tributions of the CSS, whose only non-zero density ma¬ 
trix element is poo = 1 . The fitted atom numbers N 
for different angles /3 agree within 15 % with the values 
independently measured from the shift of the cavity reso¬ 
nance. We then use the density matrix p of the heralded 
state as the only free parameter, to fit the theoretical dis¬ 
tributions gth{np) to the measured photon distributions 
g{np) along all four angles /3. We do this by minimizing 
the least squares deviation D weighted by the error <j g of 
g(n p ), given by 

9th(np,p) -g{np,p) ~ 2 

^9 

Since the photon distributions g(np) measure we 
can obtain the even terms of the density matrix (p mn 
where m + n is even) and are not sensitive to the odd 
terms. Because the overall heralding probability is pq = 
1 . 5 %, the higher-order Dicke state components are ex¬ 
ponentially suppressed. We fit the density matrix up to 
Dicke state |4) x . The fitted values p22 = 0 . 03 T 0 . 02 , p 33 = 
0.02 =b 0.01, P44 = 0.01 ± 0.01 agree with the theoretical 
expectation|l] for our system. 


^ = EE 

(3 n> 0 


From the fitted density matrix p (with coherence 
terms) we obtain the spin distributions f(Sp) using (|8| 
for different angles /?, as shown in Fig. 2e-h of the main 
text. 

To reconstruct the Wigner function for the spin state 
on the Bloch sphereJTJ [4], we convert p from the Dicke 
state basis into the spherical harmonic basis and obtain 
the normalized Wigner function according to 

N k 

EE PkqYkq{0A)i (l 2 ) 

k =0 q=—k 

where the terms pk q represent the density elements in the 
spherical harmonic basis and Yk q (0,(j)) are the spherical 
harmonics, with #, (j) being the polar and azimuthal an¬ 
gles on the Bloch sphere respectively. The normalization 
factor y/ 2S/tt is chosen such that the CSS has 0) = 
1. Note that, in the limit of large atom number, this nor¬ 
malization also means that the pure first excited Dicke 
state has 0) = —1, and generally the value of the 

Wigner function on the x axis depends only on the pop¬ 
ulations p nn such that W(0 = |, 0 = 0) = J] n (-1 ) n p nn - 

MEASUREMENT OF MEAN VALUE OF S z 

The measured photon distributions g(np) do not give 
information about the density matrix odd terms (p mn 
where m + n is odd). In order to bound the odd terms 
we verify that the heralding process does not displace the 
produced heralded state relative to the CSS. This is ac¬ 
complished by performing a measurement with a probe 
beam polarized at 45 degrees relative to |v), such that 


w(o,4>) = 


VzsJ 
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the difference between the measured |h) and |v) photon 
numbers is proportional to S z . We find a heralding- 
light-induced shift S(S Z ) = —0.2 ± 1.6, consistent with 
zero, and very small compared to the CSS rms width 
(AS^css ~ 30. Therefore we set the odd terms of the 
density matrix to zero in Fig. 3b-d. 


ENTANGLEMENT DEPTH FOR FINITE 
CONTRAST 


Y^Wi = 1. For each PN-i, 

B(P 0 ,k,N-i) = B(P 0iN -i,k,N-i)/(N-i) 

< B(P 0iN -i,k,N)/(N-i). (15) 

Here Po,at-z is the probability for the state to be found 
in the ground state in the subspace of total spin N — i. 

Measurements of the spin distributions do not allow us 
to determine the total spin of the system at single-atom 
resolution. We define populations of the CSS and the 
first Dicke state by 


Entanglement depth is defined as the minimum num¬ 
ber of entangled particles in an ensemble. A fully sepa¬ 
rable pure state can be written as \cp) = \(fi) 0 ... (8) \(Pn), 
where N is the atom number. A pure ^-producible state 
can be written as \<p) = \p>i'"' kl ) ® ••• <8) | p]^" ,kM )^ where 
fci,..., kM < k, ki + ... + kM = N. If a state cannot 
be written as a pure (k — l)-producible state or a mixed 
state of (k — l)-producible states, then it has entangle¬ 
ment depth of at least k. 

We slightly generalize the entanglement criterion de¬ 
rived in Ref. [5] to take into account the finite contrast 
C of the collective atomic spin in our experiment. The 
derivation in Ref. [5] considers the case in the fully sym¬ 
metric Dicke subspace of N atoms, and finds that for a 
^-producible state the maximum population of the first 
Dicke state pn (Pi) as a function of the CSS population 
Poo (Po) is 


max Pi 
Pq 


Po 

N 


max 


Vk max Pm-i(«i, ..., clm-i) 


n 


+ VtF 1 (y/P 0 /x) 


(13) 


Here M = [N/k ], k' = N - k(M - 1), and 

F n (ai,... ,a n ) = J2i =1 AP • Equation (13) is gener- 
ally not a concave function of Po- In order to obtain the 
upper bound for mixed states, denote the concave hull of 
the right side of equation ( [l3| ) as S(Po, k , N). We define 
B(P 0 ,k,N) = B(P 0 ,k,N)/N. Note that when Ai < N 2 , 
S(P 0 ,*,JVi) <B(P 0 ,k,N 2 ). 


The heralded state we produce does not necessarily 
retain perfect contrast, so the state can be a mixture 
of different total spins S = N, N — 1,..., N(1 — e), with 
e ^ 1%. The contrast loss is mainly caused by the deco¬ 
herence between P = 1 magnetic sublevels, and the free 
space scattering of the heralding light by the atoms. We 
decompose the density matrix p into the total spin basis 


eN 

P = P J W i p N -i- (14) 

i =0 

Here pN-i is the density matrix in the subspace of total 
spin S = N — i, Wi is the weight for each pN-i and 


eN 

Po = P^WjPo^N-j, (16) 

i =0 
eN 

Pi = Y, W i P ^N-i- (17) 

i =0 

The upper bound of Pi is given by 


eN 

max Pi < Wi max Pi w-i 
P 0 “ P0,N-i ’ 


eN 

< ^2 mB(Po tN -j, k,N — i)/(N 

i =0 


0 - ( 18 ) 


Using equation (15) and the fact that R(Po,fc,AT) is a 
concave function of Pq we have 


eN 

max Pi < Y, rejP(P 0; Ar-i, k, N)/(N - eN) 

0 i =o 

1 eN 

~ (l- e )AT E (E^ P o,iv-i,fc,lv) 

= ±B(P 0 ,k,N). (19) 

Here C is the contrast of the collective spin. Comparing 
to Ref.[5], the result is modified by a factor 1/C. In our 
experiment, C = 0.99to‘o2> so the effects of finite contrast 
on entanglement depth are minimal. 


ENTANGLEMENT DEPTH IN TERMS OF THE 
ACTUAL ATOM NUMBER 

In the experiment the atoms have spatially varying 
coupling to the probe light. However, the criterion in 
Ref. [5] is derived for the case where atoms are equally 
coupled to the light. Here we generalize the entangle¬ 
ment criterion to our experimental conditions and prove 
that the sample-averaged fractional entanglement depth 
for the ensemble containing 3100 actual non-uniformly 
coupled atoms is the same as that of 2100 uniformly cou¬ 
pled effective atoms. Consider an ensemble of N a actual 
atoms where each atom j has spin component f z j and 









cooperativity rjj. The effective total spin of the ensemble 
is S z and the effective cooperativity is 77 , so that 

N a 

S zV = J2f zJ xr 1j . (20) 

3 = 1 

As mentioned in the main paper, the ideal heralded 
state |^i) (the first Dicke state of non-uniformly coupled 
atoms) is the destructive interference of two slightly dis¬ 
placed CSSs | ±0) and can be written as 


IV’i) \4>) ~ \-<t>) 


e iS zV r/(4A) _ e -iS zV r/(4A) 


1 ^ 0 ) 


„ir/(4A) Ef=“i _ „-»r/(4A) Ef=°i fzjVj 


IV’o), 

( 21 ) 


where |^o) is the initial CSS along x. By expanding the 
exponent to first order and using f z = (/ +jX — /_ ?x )/(2i), 
we get 


( Na \ 

“ 1/2 N a 

- 

1^1)= XX 2 

XV 

1 Qj')x 

W J 

3 = 1 



( 22 ) 


where | 0 j) x and |1 j) x are the single-particle spin eigen¬ 
states along x of the atom j. For a fully separable 
state | ip) = nJL tt i + Pj |lj)x + • • •) the popula- 


tion Pi = |(</#i )| 2 is given by 

( n « \ 1 

N a 

p i= tX^ 

x^ n a i' 

\i=i J 

i =1 


The expression for Pi is similar to that in Ref.|5 and 
differs by the additional weight factor r\j. When the real 
atom number N a 1, the upper bound of Pi for the 
fully separable state \tp), S(Po, A a ), as a function of the 
population Pq = |((/?| / 0o)| 2 , is the same as Ref.[5], and 
independent of N a . 

Next consider a state which can be factorized into two 
subsets | ip) = < 8 > |where k% + & 2 = N a . 

Each |can be expanded as 

\(fi) = CLil'ipQ 1 ) + bi\‘ipi z ) + (24) 

where IV’o*) is the CSS containing ki atoms, and {ipi*) 
is given by ?? with N a replaced by ki. The populations 
P 0 = |(^#o )| 2 and Pi = |((/#i )| 2 are given by 

Po = |ai| 2 |a 2 | 2 , 


/ N a \ 


fel 

N a 

£■$ 

«2^1a 

X^ 2+a 

E 

v-> / 

\ 

i=i \ 

j=k X + l 
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The expression for Pi recovers that of Ref. 5 when 
r)j = 1. When fci,fc 2 and N a are large, we take the en- 
semble averages Vj = h(v 2 ), T,f= kl +i rf = M 7 ? 2 ) 

and Y^j=i Vj = V (jf) ■ Therefore the bound of Pi 
in equation (??), B(Po,k a = max{&q, }, N a ), is the 
same as B(Po,fc,7V) for uniformly coupled atoms when 
k a /N a = k/N. This proves that the average fractional 
entanglement depth for the ensemble containing 3100 ac¬ 
tual non-uniformly coupled atoms is the same as that of 
2100 uniformly coupled effective atoms, thus in our sys¬ 
tem a minimum of 1970 out of 2100 effective atoms or 
2910 out of 3100 real atoms are mutually entangled. 

It might seem as if the addition of N w N weakly 
coupled atoms (coupling strength r] w ) to the system 
would increase the entanglement depth without having 
physical consequences as long as N w rj ^ Nr] 2 . However 
in this case the uncertainty AA f f on the entanglement 
depth also increases, given by 

so as to be consistent with the entanglement depth N 
prior to adding the weakly coupled atoms. Atoms that 
do not change the observed spin distribution have no ef¬ 
fect on the entanglement depth. 
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